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Abstract. We address nonlocality of a class of fully inseparable three-mode 
Gaussian states generated either by bilinear three-mode Hamiltonians or by a 
sequence of bilinear two-mode Hamiltonians. Two different tests revealing strong 
nonlocality are considered, in which the dichotomic Bell operator is represented by 
displaced parity and by pseudospin operator respectively. Three-mode states are 
also considered as a conditional source of two-mode non Gaussian states, whose 
nonlocality properties are analyzed. We found that the non Gaussian character 
of the conditional states allows violation of Bell's inequalities (by parity and 
pseudospin tests) stronger than with a conventional twin-beam state. However, 
the non Gaussian character is not sufficient to reveal nonlocality through a 
dichotomized quadrature measurement strategy. 



1. Introduction 

Einstein, Podolsky and Rosen (EPR) formulated their famous argument about the 
completness of quantum mechanics in the framework of continuous variable systems 
pp. However, after that Bohm gave a dichotomized version of it the debate 
concerning nonlocality moved to systems described by discrete variables, leading Bell 
to formulate his celebrated inequalities in a dichotomized fashion 3 . Recently, the 
increasing importance of continuous variable systems leads many authors to explore 
the nonlocality issue in its original setting, where dichotomic observables to test Bell's 
inequalities are not uniquely determined. The attempts to translate Bell's inequalities 
to continuous variable systems clarified the fact that crucial in a nonlocality test is the 
existence of a set of dichotomized bounded observables used to perform the test itself, 
from which the so called "Bell operator" is derived. The more debated question has 
dealt with the nonlocality of the normalized version of the original EPR state, i.e. the 
twin-beam (TWB) state of radiation produced by spontaneuos downconversion in a 
parametric amplifier 0]. Nonlocality of the TWB state was not clear for a long time. 
Using the Wigner function approach, Bell argued that the original EPR state, and as 
a consequence the TWB too, does not exhibit nonlocality because its Wigner function 
is positive, and therefore represents a local hidden variable description [5]. More 
recently, Banaszek and Wodkiewicz 6 showed instead how to reveal nonlocality of 
the EPR state through the measurement of displaced parity operator. Furthermore a 
subsequent work of Chen et al. showed that TWB's violation of Bell's inequalities 
may achieve the maximum value admitted by quantum mechanics upon a suitable 
choice of the measured observables. Indeed, the amount of violation crucially depends 
on the kind of Bell operator adopted in the analysis, ranging from no violation to 
maximal violation for the same (entangled) quantum state. 
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Systems which involves only two parties are the simplest setting where to study 
violation of local realism in quantum mechanics. A more complex scenario arises 
if multipartite systems are considered. Studying the peculiar quantum features 
of these systems is worthwhile in view of their relevance in the development of 
quantum communication technology e.g. to manipulate and distribute information 
in a quantum communication network [HI E]. Although the study of multipartite 
nonlocality has originated without the use of inequalities an approach to 

derive Bell inequalities has been developed a l so f° r these systems and applied 
to characterize their entanglement properties |12| . Being originally developed in 
the framework of discrete variables, these multiparty Bell inequalities have found 
application also in the characterization of continuous variable systems |13l 114) . 

The aim of this paper is to apply the various approaches hitherto developed to 
test nonlocality of two- and three-mode continuous variable systems. We will consider 
tripartite Gaussian states as well as non-Gaussian bipartite states. In the first case 
strong violation of Bell inequalities is found, allowing the Bell factor to reach values of 
B ~ 3, while in the second case enhancement of nonlocality is obtained in comparison 
with the TWB case. 

The paper is organized as follows. In Sec. [3 we review the different approaches 
to test nonlocality in the framework of continuous variables and introduce notation 
that will be used throughout the paper. The three-mode states we are interested in 
are introduced in Sec. |3J and their violation of local realism is analyzed in 13. II and 13. 21 
In Sec. 01 the tripartite states are considered as sources for conditional generation of 
non-Gaussian bipartite states, whose nonlocal proprties are then studied in 14.11 14.21 
and 14.31 Finally, the main results obtained are summarized in Sec. El which closes the 
paper with some concluding remarks. 

2. Nonlocality tests for continuous variables 

In this section we will briefly recall the inequalities imposed by local realism in the 
cases of our interest. Let us start by focusing our attention on a bipartite system. 
Let m{a\) = ±1 and m(a' 1 ) = ±1 denote two possible outcomes of two possible 
measurements on the first subsystem and similarly m(a 2 ) = ±1 and m(a' 2 ) — ±1 for 
the second subsystem. The essential feature of this measurements is that they are 
local, dichotomic and bounded. The Bell's combination 

F 2 = rn(ai)m(a 2 ) + m(ai)m(a' 2 ) + m(ati)rn(a 2 ) — m(a' 1 )m(a 2 ) (1) 

under the assumption of local realism gives rise to the well known Bell-CHSH 
inequality |15| : 

B 2 = \E(a 1 ,a 2 )+E{a 1 ,a' 2 )+E{a' 1 ,a 2 )-E{a' 1> a' 2 )\<2, (2) 

where E(a±, a 2 ) is the correlation function between the measurement results, i.e., the 
expectation value of the products of the results of the experiments m(ot\) and m(a 2 ). 

In the case of a n-partite system, a nonlocality test is possible using the Bcll- 
Klyshko inequalities ^2 ^] which provides a generalization of inequality J2J • These 
inequalities are based on the following recursively defined linear combination: 

F n = - [m(a n ) + m(a' n )] F n _i + i [m(a„) - m(a' n )] F^ l _ 1 , (3) 

where m(a n ) — ±1 and m(a' n ) = ±1 refer to measurements on the n-party of the 
system, and F' n denote the same expression as F n but with all the ay and a'j inverted. 
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In the case of a three-partite system, local realism assumption imposes the following 
inequality from combination |J3J: 

B 3 = \E(ai, ct2, a' 3 ) + E(ai, a' 2 , as) + E(a' lt a 2 , a 3 ) — E(a[, a 2 , a' 3 )\ < 2 , (4) 

where again E(a\, a-z, a 3 ) is the correlation function between the measurement results. 
Quantum mechanical systems can violate inequalities J2J) and(@J) by a maximal amount 
given by, respectively, B 2 < 2^/2 and Z? 3 < 4 (see, e.g., Ref. [T2]). 

We now briefly review three different strategies to reveal quantum nonlocality in 
the framework of continuous variables systems. Recall that in the case of a discrete 
bipartite system, for example a spin-i two particle system, the local dichotomic 
bounded observable usually taken into account is the spin of the particle in a fixed 
direction, say d . Hence the correlation between two measurements performed over 
the two particles is -E(di,d 2 ) = (di<r ® d 2 er), where the operator er = (a x , a y , a z ) 
is decomposed on the Pauli matrices base and di,d 2 are two unit vectors. The Bell 
operator is then given by the expression: 

-E?2,sp = di<r ® d 2 er + d^er ® d 2 er + d^cr ® d' 2 er — d^er Cg) d' 2 cr . (5) 

Consider now a n-partite continuous variable system identified by the creation 
operator at and the annihilation operator aj (J = 1, . . . , n) with boson commutation 
relations associated. Following the original argument by EPR it is quite natural 
attempting to reveal the nonlocality of this system trying to infer quadratures of 
one subsystem from those of the others. From now on, we will refer to this procedure 
as a "Homodyne nonlocality test" , as quadrature measurements of radiation field 
are performed through homodyne detection. Here we identify the quadrature Xj (0) 
according to the definition X® = -4= (aje~ z9 + a)-e %e ). As they are local but neither 
bounded nor dichotomic, quadrature observables are not immediately suitable to 
perform a nonlocality test based on Bell's inequalities. The procedure to make 
them bounded and dichotomic is quite arbitrary and consist in the assignment of 
two domains D + and D_ to each observable JH]. When the result of a quadrature 
measurement falls in the domain D± the value ±1 is associated to it. Usually the choice 
D± — M. is considered, though a choice suitable to the system under investigation 
may be preferable. Considering a bipartite system we can introduce the following 
quantities: 
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(6) 



where P(xi,x 2 ) is the joint probability distribution of the quadratures x\ and x 2 . 
We can now identify the homodyne correlation function Eu{9, <p) as 

E H (9, V ) = P ++ (xlx^+P..(xt,xt)-P + M > x^-P. + (x 9 1 ,x^) ,(7) 

which can be straightforwardly used to construct the Bell combination B 2 .h of Eq.JSJ 
and to perform the nonlocality test. The main problem of pursuing such a nonlocality 
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test is that it is not suitable in case of systems described by a positive Wigner 
function, as the TWB state of radiation defined as \X) = y/l — X 2 J2 n X n \nn), where 
X = tanhr and r is the squeezing parameter. Indeed, a positive Wigner function can 
be interpreted as a hidden phase-space probability distribution, preventing violation 
of Bell-CHSH inequality unless the measured observables have an unbounded Wigner 
representation, which is not the case of the dichotomized quadrature measurement 
described above. Considering in fact that P(x 6 1 ,X2) can be determined as a marginal 
distribution from the Wigner function one can write from Eqs. © and J7J: 

E H (6,tp)= f dx\dxldx 6 ^ dx^ sgn{x\,xl)W(x\,x^ ,x%,x% + %) , (8) 

where the integration is performed over the whole phase-space and without loss of 
generality we have considered D± = M ± . Eq. © itself is indeed a local hidden variable 
description of the correlation function, hence obeying inequality • 

In order to overcome this obstacle different strategies have been considered by 
many authors, based essentially on parity measurements. Banaszek and Wodkiewicz 
[HI demonstrated the nonlocality of the TWB considering as local observable on 
subsystem j the parity operator on the state displaced by ay (hence we will refer 
to this procedure as a "Displaced Parity (DP) nonlocality test"), which is dichotomic 
and bounded: 

11(a) = ^DjiajX-irW^). (9) 

In the above formula, a — (a%, a n ), while nj — a) a and Dj(otj) — exp[aja}- — a*jaj] 
denote the number operator and the phase space displacement operator for the 
subsystem j. Hence the correlation function reads: 

E DP {a) = (11(a)), (10) 

from which Bell's combinations &2,dp m Eq.J2J and B^^p in Eq.Q) can be easily 
reconstructed in the cases n = 2, 3. The reason why this procedure would be able to 
reveal nonlocality also in case of quantum states characterized by a positive Wigner 
function is clear using the following relation: 

W(a) = (~Y (11(a)) . (11) 



(12) 



Indeed, the analog of Eq.lJSJ is 



E DP (a) = J d 2n X W(a)5 {2n \a - X) 



Being the Dirac-c) distribution unbounded, Ineqs. @ and (0J) are no more necessarily 
valid for $2, dp and S3, up. The maximal violation found with this procedure for a 
EPR state is B2.DP — 2.32 ^7], still far from the maximum violation admitted by 
quantum mechanics. 

Another strategy, developed by Chen et al. [7], shares a similar behavior as the 
one described above, allowing to reveal nonlocality for quantum states with positive 
Wigner function. Interestingly, this type of nonlocality test, which we will refer to as 
"Pseudospin (PS) nonlocality test" , admits a maximum violation for the EPR state. 
It can be seen as a generalization to continuous variable systems of the one introduced 
by Gisin and Peres for the case of discrete variable systems ^HJi hence, for the case of 
a pure bipartite system, it is equivalent to an entanglement test |17j . Let us consider 
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the following set of operators, known as pseudospins in view of their commutation 
relations, = (s 3 x , s 3 y ,sl) acting on the j-th subsystem 

oo 

= £(|2n+l) J -<2n + l|-|2n) J -<2n|), 

n=0 

4±4 = 24, 

d j s j = s j z cos^' + sm.6 j (e ivi sL + e~ 1 ^ s J + ), (13) 

where = |2n)j(2n + 1| = {s+Y and cP is a unit vector associated to the 

angles J and <pi , In analogy to the spin-i system and defining d = (d 1 , ...,d n ) the 
correlation function is simply given by: 

E PS (d) = (®* =1 d J V) , (14) 

from which the Bell combinations Z?2,ps and £>3,ps are evaluated. Also different 
representations of the spin-i algebra have been discussed in the recent literature 
^3 GO]- In particular in Ref. [20] it has been pointed out that different representations 
lead to different expectation values of the Bell operators. Hence, the violation of Bell 
inequality for continuous variable systems turns out to depend, besides to orientational 
parameters, also to configurational ones. In the following sections we will also consider 
the pseudospin operators IF = (IGj., 11^, IPj taken into account in Ref. [20|, which have 
the following Wigner representation: 

W ui =sgnx j w nl = -S(x 3 ) V- W n ,=-«6{xj)5(3,j) > (15) 

where Xj = Xp yj = Xj and V stands for the "principal value". The correlation 
function obtained using operators IF will be indicated as E' PS (d) = d J IF). 



3. Three-mode nonlocality 

In this section we will analyze the nonlocal properties of tripartite Gaussian states. 
In particular we will consider two classes of states, the first one proposed by Van 
Loock and Braunstein ^21, the second one proposed in |2T]. The reason why we 
consider this two classes is that the first is a very natural and scalable way to produce 
multimode entanglement using only passive optical elements and single squeezers, 
while the second one is the simplest way to produce three mode entanglement using 
a single nonlinear optical device. Indeed, both states can be achieved experimentally 
\22\ \2'.'>\ . As concern the first class of states, it is generated with the aid of three 
single mode squeezed states combined in a "fritter" (a three mode generalization of 
a beam-splitter). The evolution is then ruled by a sequence of single and two mode 
quadratic Hamiltonians. As a consequence, being generated from vacuum, the three- 
mode entangled state is Gaussian and its Wigner function is given by: 



Ws{x,y) = -^exp 



-(x,y)CT 



(16) 



where x = {x\, X2, X3), y = (yi, 2/2,2/3) are the positions and momenta of the three 
modes and C _1 is the inverse of the covariance matrix, whose explicit expression 
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(17) 



sinh 2r, S 



| cosh r sinh r and r is 



cosh 2r 

the squeezing parameter (with equal squeezing in all initial modes) . The second class 
of tripartite entangled states is generated in a single non linear crystal through the 
following interaction Hamiltonian: 



Hi, 



7i a\ 4 



,t. 



B.c. 



(18) 



Hi n t describes two interlinked bilinear interactions taking place among three modes 
of the radiation field coupled with the support of two parametric pumps. It 
can be realized in media by a suitable configuration exposed in Ref. |2.'ij . 

Notice that the same dynamics can be implemented in different physical systems, 
including optomechanical couplers and Bose-Einstein condensates in the linear regime 
|24l 1251 12(i| . The effective coupling constants jj, j — 1,2, of the two parametric 
processes are proportional to the nonlinear susceptibilities and the pump intensities. 
If we take the vacuum |0) = |0)i <£> |0) 2 <£> 1 0) 3 as the initial state, the evolved state 
|T) = e~ tHintt \0) belongs to the class of the coherent states of SU(2, 1) and it reads 
[2H EZ1: 

p/2 



1 v ( n 2 / n 3 \ 



q/2 



^-i(p<p2+q<p3) 



(P±fl|p+g,p, ff ),(19) 

plq] 



where Nj(t) = (a]j(t)a(t)) represent the average number of photons in the j-th mode 
and <pj are phase factors. The explicit expressions of Nj(t) are: 

l7i| 2 |72| 2 



N 2 = 
N 3 = 



r> 4 
I71I 2 



[cos rtt — iy 



f2 2 



sin 2 (m) , 



(20) 



with f2 = A/I72I 2 — |7i| 2 and Ni = N2 + N3. Also for this second class, being the initial 
state Gaussian and the Hamiltonian quadratic, the evolved states will be Gaussian. 
The Wigner function reads as follows |21l l2"%| : 



Wr(x,y) = -^exp 



-(x,y)V- 



(21) 



where V 1 is the inverse of the covariance matrix, whose explicit expression is: 
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where 

A = 2^2(1+^1)008 02 V = 2y/N 2 (l + N 1 ) sin 2 T = 2N 1 + 1 

6 = 2^3(1 + ^1) cos 03 S = 2^3(1 + ^1) sin 3 = 2iV 2 + 1 

C = 2v/iV^V3COs(0 2 -0 3 ) £ = 2v/iV2~JV3sin(02-0 3 ) = 2A^ 3 + 1 . (23) 

Both classes of states are fully inseparable for any value of the coupling constants, 
namely cannot be written as a factorized state for any grouping of the modes. 
Therefore they are good candidates to reveal true tripartite nonlocality. Being 
Gaussian states, however, nonlocality cannot be revealed by homodyne detection. 
In the following we analyze the results for displaced parity and pseudospin tests. 

3.1. Displaced parity test 

Let us start the study of tripartite systems nonlocality using the "displaced parity 
test". Considering the correlation function Ejjp{ol) given by Eq. i|l(J|) . the state 
l |lti| l was found in |13j to give a maximal violation of B 3- dp — 2.32 in the limit of 
large squeezing and small displacement. The study in |13| however was performed 
for a particular choice of displacement parameters: a% = a 2 = a 3 — and 
a 'i — a 2 = a 3 — i^f~3 ■ A numerical optimization of the displacement parameters lead 
us to identify a number of parameterizations that allow a significantly higher violation 
of Bell's inequality. As an example, consider the one given by a.\ = a 2 = ct 3 = i\fj 
and a'j = a' 2 = a' 3 = from which follows that 

B 3 , DP = 3cxp (-12e- 2r j) - exp (24e 2r j) , (24) 

hence the remarkably high asymptotic value of B^^dp — 3 is found for large r and 
J ^ (see Fig. . The importance of a suitable choice of the displacement parameters 
is apparent if this asymptotic value is compared to the violations obtained in the 
nonlocality study performed in Ref. 13 . In that work in fact generalizations to more 
than three modes of state (|16f) were also considered, giving an increasing violation 
of Bell inequality as the number of modes increases, but never founding a violation 
greater than 2.8 . Determining the optimal choice of the displacement parameters for 
a given state is in general a challenging task. To our knowledge indeed there exists 
no general prescription to find it out, and ultimately one must rely onto a numerical 
analysis (see, e.g. EH])- Nevertheless, a careful inspection of the symmetries of the 
state under consideration may be helpful. In order to clarify this observation let us 
consider the explicit expression of the correlation function Ep,p{ct\, a 2 , 013) for state 

as, a 3 ) = exp j - ^e 2r [(yt +y 2 + y^f + (x 2 - x 3 ) 2 + (x 2 - xi) 2 + (xi - x 3 ) 2 ] 

~\ e ~ 2r + ^2 + x 3 ) 2 + (y 2 - y 3 ) 2 + (y 2 - Vl ) 2 + ( yi - y 3 ) 2 ] 1(25) 

Hence, from Eq. @ it follows that the Bell combination B 3 _dp is given by the sum 
of three positive and one negative term. It is reasonable to expect that the maximal 
violation of nonlocality will be achieved for large r. We see from Eq. Ij25(l that, in 
this limit, all the correlation functions in B 3 ,dp vanish for nonzero displacements, 
unless the coefficients of e 2r are zero. Hence we impose the following system of 
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equations, which allows for the three positive terms in B^jyp not to vanish (we consider 
ot-k = x k +iyk and a' k = x' k + iy' k for A; = 1, 2, 3): 

(Vi + V2 + Vs) 2 + Oa - x 3 ) 2 + {x 2 - x\) 2 + (x[ - x 3 ) 2 = 

(m + y*2 + V3) 2 + (x' 2 - x 3 ) 2 + {x' 2 - xi) 2 + (xi - x 3 ) 2 = 

(2/1 + 2/2 + y' 3 ) 2 + {X2 - 4) 2 + (*2 - x x ) 2 + ( Xl ~ x' 3 ) 2 = . (26) 

We see that the parameterization used to obtain Eq. (|24f> is a solution of this system. 
Clearly, any other solution will give, in the limit of large r, the same violation given 
by Eq. (|24|) . namely B 3 dp — » 3. In order to compare the violation of Bell's inequality 
admitted by the state with the one that will be obtained below considering the 
state l|19[) j it is useful to rewrite Eq. I|24l) as a function of total mean photon number 
N = Ni + N2 + N 3 . Given that N = 3 sinh 2 r and optimizing the displacement J 
we obtain the result shown in Fig [3] The asymptotic expression of the optimized 

displacement as a function of N is J = g^-ArcSinh [y^Pj, hence very small angles 
are required. 




Figure 1. Plot of the Bell combination 1241 . Only values violating Bell Inequality 
fll are shown. 



Consider now the tripartite state generated by the Hamiltonian l|18fl . The 
correlation function is now given by Eq. I|1U|) through the Wigner function l|21|) . 
The symmetry of the state suggests a maximum violation of Bell inequality for 
N2 = N 3 = (recalling Eq. I|20|ll. while the fact that the separability of the state 
doesn't depend on the phases </>2 and <fo [2] suggests that they are not crucial for the 
nonlocality test. If we consider the same parametrization that led us to Eq. I|24|l and 
fix cf>2 — (f>3 — 7T, we find: 

_l + e 6J ( 1 + N + 2 V2 y/N (2+JV)) J {±+7N+6s/2^/n (2+JV)) 

B3, dp = 7 , =r , (27) 

g 4 J (3+3 N+2 ^2 y/N(2+N)) 

from which follows an asymptotic violation of Bell's inequalities of B^ : dp — 2.89, for 
large N and small J . A slightly better result is found if a parametrization, more 
suitable and numerically optimized for state l|19|) . is considered: ct\ = % > \fj 1 cx2 = 
«3 = ol x = 0, a' 2 — —\f~3-, a' 3 = \fj <t>2 — and <\>2 — it. The Bell combination B% t DP 
for this choice of parameters is depicted in Fig. [21 We note that in this case a larger 
choice of angles allows the violation of Bell inequality if compared with Fig. ^ As 
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Figure 2. Bell combination obtained choosing optimized displacement 
parameters for state 1191 (see text for details). Only values violating Bell 
Inequality J1J are shown. 




Figure 3. Bell combination By. dp for state 1161 . dotted line, and 1191 . solid 
line. The displacement parameters J have been optimized to give the maximum 
value of B for N fixed. 



before, optimizing the displacement J for each N, it is possible to find the maximum 
violation of Bell inequality as a function of N. We find that the asymptotic relation 
between the optimized displacement and the total photon number is now JN ~ 3.21, 
confirming that not too small displacements are required. The asymptotic violation 
of Bell's inequality is now B-^^dp — 2.99. To compare the results obtained for the 
two states ijTT))) and <jl"9"|> we have plotted in Fig. |3]the Bell combination Bz.dp as a 
function of the mean total energy TV, while the displacement J has been chosen in 
order to maximize B^.dp at fixed energy. Notice that even if the two states have quite 
the same asymptotic violation, state 111 (it reaches it for lower energies. 

3.2. Pseudospin test 

Consider now the pseudospin nonlocality test. Let us calculate the expectation value of 
the correlation function l|14fl for the state |T) (for simplicity we consider 02 = </>3 = 0). 
The only non vanishing contributes are given by: 

ci = (si ®s 2 x ® si) = (si ® si ® si) 
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2{1 + N 1 fn 2^y i + Nl ) \i + N J (2s)!(2i)!V 2s + 1 ' 1 j 

and by (s* 65 s 2 65 s 3 ) = 1. The correlation function then, according to Eqs. (|13f) and 
l|14fl . reads as follows: 

Eps(d) — cos 9 1 cos 6* 2 cos # 3 + ci cos 6* 1 sin# 2 sin 9 3 (cos (£ 2 cos <p 3 + sin^ 2 sin ip 3 ) 
+ C2 cos9 2 sinO 1 sin 9 3 (cos ip 1 cos ip 3 — sint^ 1 sin</? 3 ) 
+ C3 cos 9 3 sin 9 1 sin # 2 (cos ip 1 cos <y9 2 + sin tp 1 sin <p 2 ) (29) 

Hence, without loss of generality, we can fix for example ip 1 — and tp 2 — p 3 = rr 
and look for the maximum violation of Bell inequality constructed from Eq. i|29|) . 
Notice that if the coefficients c 2 ; (i — 1,2,3) were equal to one then the maximum 
violation admitted, B^^ps = 4, should be reached. Considering Eqs. (|28|l two limiting 
cases can be studied. First, for large N2 and small -/V3 (or vice-versa) a numerical 
evaluation of the coefficients cj shows that C3 — > 1, while the other two vanish. Hence, 
considering 9 3 — 0, the correlation function l|29|l reduces to that of a TWB subjected 
to a pseudospin nonlocality test (see Eq. (|39|l below), hence allowing an asymptotic 
violation of B^^ps = 2V2. This result should be expected, since in this limiting case 
the state l|19JI reduces to a TWB for modes a\ and a 2 , while mode 03 remains in the 
vacuum state and factors out. Consider now the case in which N% = N3 = A 
numerical evaluation shows that the coefficients Cj — ► | for large N, hence also in this 
case the maximum violation cannot be attained. The asymptotic violation turns out 
to be Bs^ps — 2.63. 

As already mentioned in Sec.|21othcr representations for the pseudospin operators 
can be considered. Using Eqs. (fl5|l and (|2"T|) it is possible to calculate the correlation 
function E' PS (d). Setting again the azimuthal angles ip 1 = 0, the latter shows the 
same structure as Eps(d) where now the coefficient q are replaced by: 



2arctan (^F) , , 2arctanVAT 



^ n(l + N) - 2 --3- 1T(;i+ N ) ■ ( 3 °) 

An appropriate choice of angles leads to a maximal violation of Bell's inequality given 
by Bz y ps — 2.22 (see Fig. 0}, which is now reached for N ~ 1, value for which 
the coefficients c\ are approximately near their maxima. As already pointed out, we 
see that different representations of the pseudospin operators give rise to different 
expectation values for the Bell operator. 

Applying now the same procedure to state (|16f) we find the same structure for 
the correlation function E' PS , where the coefficients are now given by: 

_ 6 arctan( 4co f (, ' ) sinh(r) ) 

J - J - J - V 3 ( 2 + e4r ) ( 31 ) 



7T a/5 + 4 cosh(4r) 
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After an optimization of the angles Q l we obtain a maximal violation of B^^ps — 2.09 
(see Fig.0J, for r ~ 0.42 (N ~ 0.56) that maximizes the coefficients Cj. 



2 . 2 
2 . 1 

m 2 



1.9 
1.8 

2 4 6 8 

N 

Figure 4. Plot of Bell combination Bg^pg for states 1211 and 1161 . solid and 
dashed line respectively. N is the total number of photons as in Fig. [3] 



4. Degaussified state and two mode nonlocality 

In this section we consider the tripartite state l|19fl as a source of two-mode states. In 
particular, we study the nonlocality of a two-mode non-Gaussian state obtained by a 
conditional measurement performed on state l|19l) . Gaussian states are at the heart 
of quantum information processing with continuous variables. The reason for this is 
that Gaussian is the character of the vacuum state of quantum electrodynamics. This 
observation, in combination with the fact that the quantum evolutions achievable 
with current technology are described by Hamiltonian operators at most bilinear 
in the quantum fields, accounts for the fact that the states commonly produced in 
laboratories are Gaussian. In fact, bilinear evolutions preserve the Gaussian character 
of the vacuum state. In addition, it is worth noticing that the operation of tracing 
out a mode from a multipartite Gaussian state preserves the Gaussian character too, 
and the same observation is valid when the evolution of a state in a standard noisy 
channel is considered. Indeed, the only feasible way to "degaussify" a state is through 
a conditional measurement, or by statistically mixing it with another Gaussian state. 
The reason to study non-Gaussian states is that when the Gaussian character is lost, 
then immediately the Wigner function of the state becomes negative, for pure states, 
hence stronger nonclassical properties should emerge. Actually, various authors have 
recently investigated the nonlocality properties of two-mode non-Gaussian states. In 
particular, a twin-beam state subjected to inconclusive photon subtraction (IPS state) 
has been considered in Refs. [3D] and while in Ref. J7| it has been pointed out 
that if the entangled coherent states [32] could be produced experimentally they would 
allow for the maximal violation (i.e., B2 — 2^2) both in case of a DP test as well as 
a PS test. 

The most natural way to obtain a non-Gaussian state from a Gaussian one is by 
elimination of its vacuum component. In fact, such a state is necessarily described 
by a negative Wigner function (in fact (0|f>|0) oc J d 2n aW(a)e~ 2 ^ ). Do to the 
structure of state 1)19(1 its vacuum component can be subtracted by a conditional 
measurement on mode CI3, the same observation being valid for mode a-i- Consider 
a photodetector able to distinguish only the presence or the lack of photons, i.e., an 
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ON/OFF photo-detector, and the state Q\ conditioned to the presence of at least 
one photon. The probability operator measure (POVM) is two-valued {IIojIIi}, 
n + III = I: with the element associated to the "no photons" result given by 

no = Ii<g>l2<8^(l-77)>)33<n|, ( 32 ) 

n 

where rj is the efficiency of the photodetector. The probability of the outcome is given 

by 

Pi = Tri23 [|T)(T| IIi] = TJ^^VgJ > 
while the conditional output state reads as follows 
ei = ^Tr 3 [|T)(T|ni] 



(33) 



1 + ^3 

(l + JViJjjJVa 



St 



p=i 



iVi 



pi 



No 



1 + iVi 



n+n 



\nn) (n' n'\a p 



1 + vN 3 



(l + iVi + 7V 2 )r/iV3^ \1 + 



Ni 



fl - v) p i 
pi 



(34) 



where we have identified the TWB with X = 



N 2 



i+N ■ ^° calculate the Wigner 
function associated with state p%, consider that the characteristic function of the 
POVM IIi reads as follows: 



XpiK/x) = tt<5 2 (m) - -exp 



V 



2 V 



IMI 



hence the characteristic function of q\ is given by: 

x[ei](A 1 ,A 2 ) = i(x[|T)(T|](Ai,A 2 ,0)-i f ^ x [|T)<T|](Ai, A 2 ,/*)exp 

After some algebra the Wigner function associated with state p\ can now be calculated. 
It reads as follows: 

2 



(35) 

|2 2_V 
2?7 



Wi(x,y) 




1 



VdetT^ 



exp 



-(x,y) (V) -1 



l 



V \ir J y/detD 



exp 



■(x,y) (D- 1 )' 




(37) 



where, from now on, x = (xi, x 2 ), y = (2/1,2/2), and D = V + diag(0, 0, -^,0,0, -^). 
In order to simplify the notation we have indicated with O' the 4x4 matrix obtained 
from the 6x6 matrix O deleting the elements corresponding to the third mode (3-th 
row/column and 6-th row/column), due to the trace over the 3-th mode. Of course, 
the easiest way to obtain a bipartite state from state (|19f) is to discard a mode, say 
the third one, by tracing over it. The state qt t then obtained is simply given by the 
following Wigner function: 

2 



W T r(x,y) 



: exp 



-(x,y) (v'r 1 



(38) 



,71-/ 4Vdet V 

Being the Wigner function Wtt Gaussian, we expect that this state will exhibit weaker 
nonlocality with respect to state i|3~4"|) . In the rest of the section the nonlocal properties 
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of the usual TWB state and of the states (|34[1 and l|38|) will be compared. Notice that 
state (|16|1 can be considered as an extension to three modes of the TWB. All the three 
nonlocality tests introduced in Sec. [3 will be taken into account. 

4-1. Displaced parity test 

We first study the violation of inequality (J2J) in the case of a "displaced parity 

test". As already mentioned, in Ref. [H| Banaszek and Wodkiewicz found for the 

first time that the TWB state exhibit a violation of local realism. They obtained the 

following asymptotic violation for infinite energy: Bi,dp — 2.19. Generalizing their 

procedure this result can be improved, yielding to a maximum asymptotic violation 

of B2.DP — 2.32 |17|. We have considered the following parametrization to obtain the 

maximum violation for a TWB: ct\ = —a 2 — i^fj and a[ = —a' 2 — —2>i\[3 '. The 

asymptotic relation between the squeezing parameter and the displacement angles is 
e 2rj = lo|3 Uging the 

same parametrization and considering the Bell combination 
$2,d p for the state QTr, it turns out that the same asymptotic value of the TWB 
is reached for large N 2 and small N3. In fact, as already noticed, when this limit 
is considered the original tripartite state (|19|) reduces to a factorized state composed 
by a TWB and the vacuum state. Consider now the conditional state p\ and again 
the case of large N 2 and small N3, say N% = lCP 2 -^-. As in the tripartite case 
the phase coefficients fa and fa play no rule in the characterization of nonlocality. 
A stronger violation of Bell inequality is then found and it is depicted in Fig. JSJ), 
where the parametrization ct\ = \ct 2 = \ a 'i = i\f7f & n d ct' 2 = has been adopted. 
Indeed the asymptotic violation is higher then the previous, namely B2DP — 2.41. 




Figure 5. Bell combination obtained choosing optimized displacement 
parameters for state q\ (see text for details). Only values violating inequality 
Ol are shown. 

It can be found, for large N 2 , when JN 2 ~ 0.042. A comparison with the violation 
of nonlocality attained with a IPS state shows an identical asymptotic behavior |31| . 
Nevertheless, the scheme proposed here offers the advantage that the production rate 
of state Qi [i.e., the conditional probability l|33l) ] is much greater then the production 
rate of IPS state [see Ref. 123) Eq.(14)]. This is due to the fact that only a single 
ON/OFF detection is required to produce g±, rather than the coincidence of two 
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ON/OFF detections for the case of IPS state. This could be useful from a practical 
viewpoint. 

4-2. Pseudospin test 

Let us now focus on the "pseudospin nonlocality test" . Considering a TWB state, it 
is known that the correlation function l|14(l has the following expression (setting to 
zero the azimuthal angles) [7]: 

Eps (01,62) = cos 61 cos 6 2 + Jtwb sin 81 sin 9 2 , (39) 
where, denoting with N the total photon number, 



Jtwb 



(40) 



1 + N 

It turns out that the violation of Bell inequality in this contest increases monotonically 
to the maximum value of 2y2 as the function Jtwb goes to unity. A straightforward 
calculation shows that an expression identical in form to Eq. (|39Jl can be found both 
in case of state q\ and QTr, where the following functions /1 and fr r can be identified: 



/i = 2 



far 



N 2 (1 + iV 3 r)) 



1 + Ni N 3 (1 + Ni) r) 



N 2 



1 



1 + JVi 1 + N x ^ n 

p,q=0 



E 

k,p=0 
1 + N! 



(2k+p) \ j2k + p+l 
(2k)\p\ 



2p 



N, 



2k + 1 

2q 



No 



1 + JVi 



(2p 



(1-(1- V ) P ) 
+ 2q)\ [2q~+2p+l 



1 + NJ (2p)!(2g)! 



2p+l 



In order to compare the violations in the three different cases, let us fix as in the 
previous subsection a small value for N 3 . A plot of the functions Jtwb, fi and fx r 
versus the total number of photons of the TWB for the former and of the initial three- 
partite state for the latter two is given in Fig. © . It can be seen that state Q\ achives 
large violations for smaller energies with respect to the other two states. Finally, a 
comparison with the violation attained with the IPS state may be found in Ref. |84| . 
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Figure 6. Comparison between the values of the functions /twb (solid line), fl 
(dotted line) and fx T (dashed line) defined in the text (the summation has been 
numerically performed for rj = 0.8 and ./V3 = 0.1). 
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4-3. Homodyne detection 

The negativity of the Wigner function (|37H may suggest to perform a nonlocality 
test based upon a Homodyne detection scheme. While the positivity of a Wigner 
function avoid to violate Bell inequality with such a test, its negativity is yet 
not sufficient in general to ensure a violation. Quantum states with negative Wigner 
function that doesn't violate local realism with a homodyne test are given for example 
in Refs. Considering state l(34jl it is necessary to calculate the correlation function 
|JS}. Substituting the Wigner function 1|37|) into Eq. J8J and performing the integral 
we obtain the following result: 



EhW 



1 + yNa 
4r?7V 3 



7T/ VdetZ? 



VdetV' 

2ir (-I + N3 (-2 + 7?)) 
1 + N 3 n 



2 (1 + 2 N 3 ) 7T arctan 



2 cosi/> 



/ (l+2iV 1 )(l+2jV 2 ) _ 4 2 , 
(1 + JVi) N 2 * Luh V , 



arctan 



2 cost/" 



I (1+2 N1 — N3 rj) (1+2 jV 2 +JV3 ??) 



4 cos 2 ip 



where tp = 9 + if + <j) 2 . A plot of the correlation function l|4*2t is depicted in 
Fig. 14.3(1 for unitary efficiency 77, together with the classical correlation function of two 
spin-i particles (see caption for details). Unfortunately, the comparison shows 
clearly that the correlation given by Eq. I|42|) is always lower then the classical one, 
hence despite the negativity of Wigner function i|37|) no violation of Bell inequality is 
achievable with this scheme. 




Figure 7. Comparison between the correlation functions obtained from two spin- 
2 particles classically correlated (solid line) and from Eq. 1421 : N2 =0.5 (dotted 



line),Af 2 = 1 (dot-dashed line), 
-/V3 = 0.5 and r) = I. 



5 (dashed line). In all cases we have fixed 



5. Conclusions 

A detailed analysis of the nonlocality properties of multipartite continuous variables 
systems obtained by parametric optical systems has been presented, using the more 
recent approaches developed to this aim. We have considered in particular two 
classes of tripartite Gaussian states that seems promising for quantum communication 
purposes in order to implement multipartite quantum protocols. The results show 
that for these states a nonlocality test based on displaced parity measurements is 
more suitable to reveal violation of local realism than one based on pseudospin 
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operators. This results are just the opposite of what have been obtained for the 
bipartite case. Notice, however, that a systematic approach to pseudospin operators 
for continuous variables haven't been developed yet, hence we have only used the 
two inequivalent configurational parameterizations more suitable for calculations. For 
displaced parity test we obtained a remarkably high asymptotic value for the Bell 
parameter, Bz,dp — 3. In this case the choice of a proper parametrization, suitable 
for the state under investigation, have revealed crucial. 

We have also explored the possibility to enhance nonlocality in bipartite systems 
considering states endowed with a nonpositive Wigner function. We investigated a 
method to conditionally produce such a state from the tripartite systems considered 
above. As expected, the Bell parameter reaches a value higher than for a TWB, namely 
&2,dp — 2.41. Also in the case of a pseudospin test an enhancement of nonlocality has 
been demonstrated, while a violation of local realism using a dichotomic quadrature 
measurement cannot be achieved. 
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